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The Mellin moments of deep inelastic structure functions at two loops 
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We perform the analytic calculation of the Mellin moments of the structure functions F2 and Fl in perturbative 
QCD up to second order corrections and in leading twist approximation. We calculate the 2-loop contributions 
to the anomalous dimensions of the singlet and non-singlet operator matrix elements and the 2-loop coefficient 
functions of F2 and Fl- Our results are in agreement with earlier calculations in the literature. 



1. INTRODUCTION 

Deep inelastic lepton-hadron scattering is one 
of the best studied reactions today. It provides 
unique information about the structure of the 
hadrons and tests one of the most important pre- 
dictions of perturbative QCD. the scale evolution 
of the structure functions ||l|,Hj3|. At the same 
time, the ever increasing accuracy of deep inelas- 
tic scattering (DIS) experiments demands more 
accurate theoretical predictions and despite all 
achievements [§J|,§,01, no complete next-to- next- 
to-leading order (NNLO) analysis for DIS reac- 
tions is available. 

The necessary perturbative QCD predictions 
at 3-loops for the anomalous dimensions of the 
structure functions and the structure function Fl 
entering in the ratio R = aL/crr, are still un- 
known, except for a number of fixed Mellin mo- 
ments of F2 and Fl |^,^. These have already 
been used in NNLO analyses |^,|lO|jl^. In gen- 
eral however, this limited information about some 
fixed Mellin moments is not sufficient to allow for 
NNLO analyses of all data of DIS experiments 
and as a consequence, there are still considerable 
uncertainties on the parton densities, e.g. on the 
gluon density at small x. 

Progress beyond the state of the art has to ex- 
plore different directions, since a straightforward 
extension of rcf . to calculate more fixed Mellin 
moments, is not feasible. One such possibility 
that actually dates back to the origins of QCD 
jl}^, is to calculate the Mellin moments of the 
structure functions analytically as a general func- 
tion of N. This approach, further pioneered by 



Kazakov and Kotikov ji] to obtain the structure 
function Fl at 2-loops, turns out to be flexible 
and very promising in view of the ultimate goal, 
the anomalous dimensions at 3-loops. As a first 
step and to demonstrate the power of the method, 
we report here on the recalculation ||l^ of the per- 
turbative QCD contributions up to 2-loops to the 
structure functions F2 and Fl- 

2. FORMALISM 

The optical theorem relates the DIS structure 
functions to the forward scattering amplitude of 
photon- nucleon scattering, T^i,, which has a time- 
ordered product of two local electromagnetic cur- 
rents, jfj.(x) and jv{z), Fourier transformed into 
momentum space, to which standard perturba- 
tion theory applies. The operator product ex- 
pansion (OPE) allows to expand this current pro- 
duct around the lightcone {x ~ zY ^ Q into a se- 
ries of local composite operators of leading twist 
and spin N . The anomalous dimensions of matrix 
elements of these operators govern the scale evo- 
lution of the structure functions, while the co- 
efficient functions multiplying these matrix ele- 
ments are calculable order by order in perturba- 
tive QCD. 

In this way the Mellin moments of DIS 
structure functions can naturally be written in 
the parameters of the OPE, 



l + (-l)^ 
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where z = 2,L and all even Mellin moments are 



fixed. Anuci denote the spin-averaged hadronic 
matrix elements and Ci the coefficient functions 
and the sum extends over the flavour non-singlet 
and singlet quark and gluon contributions. 

The OPE is an operator statement and there- 
fore independent of a particular hadronic matrix 
element, so that it is standard to calculate par- 
tonic structure functions with external quarks 
and gluons in infrared regulated perturbation the- 
ory. In practice, this procedure reduces to the 
task of calculating the A^-th moment of all 4-point 
diagrams that contribute to T^^ at a given order 
in perturbation theory. This will be shown in 
more detail below. 

3. METHOD 

To illustrate the method of calculating Mellin 
moments of DIS structure functions, consider the 
following diagram, which is the most complex 
topology at 2-loops, since three propagators with 
P-momentum are present. 

A b 
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{plY{{P+Pi?) 



(2) 



{pif{pimp-P3?r{piY{{p+p^?r' 

The fat lines denote P-dependent propagators 
with P^ = 0, while the Q-momentum flows from 
right to left through the diagram. Usually, there 
are only single powers of the propagators, a = 
6 = ... = 1 in eq.(|). 

We wish to compute the 7V-th moment of this 
diagram. In a naive approach one would ex- 
pand all P-dependent propagators into sums over 
P-Pi/pj using P^ — 0. Then scaling arguments 
require the final answer for the A'^-th moment to 
be proportional to [P-Q/Q"^)^ and one is left to 
calculate 2-point functions with symbolic powers 
of scalar products in the numerator and denomi- 
nator. However, this procedure leads to multiple 
nested sums, which in general are very difficult to 
evaluate. 

More sophisticated ideas have been used in 
ref. Q . There it is shown how to determine reduc- 
tion identities for a given diagram in dimensional 
regularization in such a way, that is possible to 





set up recursion relations in the number of mo- 
ments N . A systematic classification of all 2-loop 
topologies reveals two basic building blocks. The 
two diagrams, each with only one P-dependent 
propagator. 



(3) 



are shown in diagrammatic notation similar to 
eq. (|). These diagrams can be calculated at the 
cost of one sum over F- functions. In dimensional 
regularization, Z? = 4 — 2e, the F-functions can 
be expanded in e and the sum can be solved to 
any order in e. 

The other more complex topologies like eq.(H) 
can be expressed in sums over these basic topolo- 
gies by means of recursions, so that a complex 
diagram is broken up systematically into its build- 
ing blocks. For example, the simple and beautiful 
recursion relation for eq.(y) reads 
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which relates the iV-th moment of this diagram 
to a diagram of a simpler structure, where the 
P-dependent propagator (P—ps)^ has been elimi- 
nated. The recursion eq.(Q) can be solved and for 
the second diagram on the right hand side similar 
relations are derived. 

Amazingly enough, if the recursion relations 
are set up in this way, all the multiple nested 
sums can all be solved successively in terms of 
harmonic sums, the basic functions of weight m 
being defined as follows 



(-ir 



(5) 



while higher functions can be defined recursively 

N ^ 

-'1712 , . . . ,mfc 

(0, (6) 
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Of course, this procedure requires great care 
in the way the reduction identies are applied. It 
also relies crucially on all algebraic relations for 
harmonic sums, which allow to solve the nested 
sums algorithmically and to express the result 
in the basis of harmonic sums eqs.(|^)-(Q). 

4. RESULTS 

The analytic calculation has been done with the 
symbolic manipulation program FORM ||l5|. All 
recursion relations have been implemented in a 
program, that reduces diagrams for DIS structure 
functions up to 2-loops to multiple nested sums 
over the basic building blocks and subsequently 
calls the SUMMER algorithm, to solve these 
nested sums in terms of the basis of harmonic 
sums. The database of diagrams, identical to the 
one used in ref. [|j , was generated with the help of 
QGRAF The program has been optimized 

by tabulating basic building blocks, so that the 
calculation of all 425 diagrams for F2 and Fl up 
to 2-loops is done in a few of hours on a Pentium 
Pro PC. 

The calculation is performed in dimensional 
regularization, Z? = 4 — 2e, and renormalization 
and mass factorization proceeds as described in 
refs.[|||^,^,||. To extract the gluon anomalous di- 
mension at 2-loops, we also calculate the unphy- 
sical structure functions of a scalar particle 0, 
that couples to gluons only §. These provide 
us with the necessary renormalization constants 
of the singlet operators in eq.(|]) and can be ob- 
tained from an additional operator (pF'^ F^^^ in 
the Lagrangian. 

Finally, we obtain the even Mellin moments 
of the complete set of anomalous dimensions, 
the flavour non-singlet and singlet quark and the 
gluon coefhcent functions up to 2-loops. Our 
results for the even moments of the anoma- 
lous dimensions agree with the ones published in 



refs.yj^], while our results for the coefficient func- 
tions agree with those of refs.|||,||j^J|,^. 

We have also calculated all even moments of the 
0(e)-terms of all structure functions at 2-loops, 
needed to extract the coefficient functions at 3- 
loops after mass factorization. In the appendix, 
we present our result for the even moments of the 
pure singlet quark coefficent function C2^, which 
has not been shown in this form before. The com- 
plete results and more details of the calculation 
will be published elsewhere . 
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APPENDIX 

The even Mellin moments of the pure singlet 
quark coefficient function at 2-loops. 



UfCp [ 
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130/275i(iV-l) 



-133/81 S{N^2) 
-344/275i(A^-2) 
-1714/27S'i(A^+l) + 448/275i(A^+2) 
+580/9Si{N) - 16/35i._2(A^-2) 
-f64/35i,_2(A^-l) + 64/35i,_2(A^+l) 
-16/35i,_2(A^+2) - 325i,_2(iV) 
-l-104/95i4(A^-2) - 416/95i4(7V-l) 
-272/95i,i(A^4-l) + 32/95i4(A^-l-2) 
+18A/5Si^i{N) - 16/35i,i,i(iV-2) 
+4/35i4,i(7V-l)+4/35i,i4(7V + l) 
-16/35i4,i(iV + 2) + 8S'i,i4(iV) 
+16/35i,2(A^-2) - 4/35i^2(iV-l) 
-4/35i^2(iV + l) + 16/3S'i,2(iV + 2) 
-85i,2(iV) + 5652(A^-1) 
-hl36/952(iV + l) -I- 128/9S'2(A^4-2) 
-256/352(iV) - 16524(A^-I-1) 
+16^24(^ + 2) + 852,i4(iV-l) 
-852^i,i(A^+l)-852,2(A^-l) 
+852^2(iV + l) + 253(A^-1) 
+154/353(iV + l) - 64/3S'3(A^-l-2) 



-32S3{N) - 1653,1 (A^-1) 
+ 1653,i(iV+l) + 2054(VV-l) 
-2054(A^+1)}] . 
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